Dimensionality induced entanglement in macroscopic dimer systems 
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We investigate entanglement properties of mixtures of short-range spin-s dimer coverings in lat- 
tices of arbitrary topology and dimension. We show that in one spacial dimension nearest neighbour 
entanglement exists for any spin s. Surprisingly, in higher spatial dimensions there is a threshold 
value of spin s below which the nearest neighbour entanglement disappears. The traditional "clas- 
sical" limit of large spin value corresponds to the highest nearest neighbour entanglement that we 
quantify using the negativity. 
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Consider a macroscopic system of N subsystems with 
an arbitrary spin s. Maximal entanglement of such a 
system corresponds to the configuration where pairs of 
subsytems exist in maximally entangled states (dimers). 
The amount of entanglement quantified by the relative 
entropy in this case is -jlogl^s + l)- The dimer con- 
figuration is also the most robust to noise in the sense 
that one has to destroy entanglement between each pair 
of dimers in order to destroy all entanglement. However, 
there are many dimers configuration leading to the same 
maximal amount of entanglement. Frequently in nature 
there is no reason why some of them should be more likely 
than others, which is why they occur as superpositions 
or mixtures of all possible dimer coverings. 

Here we analyze how superposing and mixing of dimers 
affects the robustness of their entanglement. We limit 
ourselves to short range dimers because in practice forces 
generating entanglement between subsystem are them- 
selves specially short range (in momentum space this, of 
course, need not be the case, a issue that is briefiy ad- 
dressed at the end of the paper). 

Recently, there has been a considerable interest in a 
possible link between entanglement and properties of 
many-body systems [H . Although a role of entanglement 
in description of such large macroscopic systems is not 
entirely clear there are some promising results showing 
the connection between the critical phenomena and en- 
tanglement [2] . Our present paper is a small contribution 
to this field. 

Let us first consider a one dimensional lattice that is a 
union of two sub-lattices La and Lb- The sites belonging 
to the lattice La {Lb) are enumerated by odd (even) 
numbers. We assume that the lattice is described by the 
so-called spin liquid, i.e., a state without magnetic order 
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with |(a6)„(„+i)) ^ -j^J2k=o"/s'\>^)n\k + a)n+i, S = 
2s + 1, a,b = 0, 1, . . . , S" - 1, 7 = exp (2hl) and the 



index n refers to the n-th site on the ring. The states 
|(a6)„(„_|_i)) form a generalized Bell basis for two parti- 
cles with spin s, i.e., they are complete and maximally 
entangled. In analogy with the spin i case a maximally 
entangled state between two sites is called a dimer and 
the states |ci),|c2) are called dimer coverings. We are 
interested in entanglement properties of the state \tp) in 
the thermodynamicl limit, i.e., for A'^ —^ 00. Note that it 
makes no difference in this limit whether we are super- 
posing or mixing two dimer coverings. As will be proven 
below all the phase information is absent from the local 
entanglement properties, which we are interested in. 

We can prove the following facts for one dimensional 
case 

1. There is always nearest neighbour entanglement for 
arbitrary value of s. Furthermore we present the 
negativity as the function of s and show that it 
increases with s. 

2. The subset of even (or odd) sites does not contain 
entanglement. 

3. The subset of even (or odd) sites is maximally en- 
tangled to the rest. 

Proof of 1. Let us first derive a density matrix for the 
nearest neighbours. Due to the translational invariance 
of the state it is enough to consider the density matrix 
P12 of the first two spins. We have 
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After some tedious but straightforward algebra we get 
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Thus, the un-normahzed density matrix pi2 reads 
Pi2^^ + \iab)){{ab)\ 

abN 

+ ^Kab)){[{a + 2Na)b]\ + h.c. (4) 

The trace of the matrix pi2 equals M — 2 + 
(,pg ^2ir^jv^ if the number of the sites is a mul- 
tipUcity of S, i.e., 2N = mS (m is an integer) and 2 
otherwise. 

In the thermodynamic hmit the normahzed state pi2 
becomes an equal mixture of the maximally entangled 
state |(a6)) and the white noise 

Pi2^'^^l^ + \\W){{ab)\. (5) 

In the similar way we can compute the density matrix 
between next nearest neighbours, for instance, between 
the first and the third spin 

Pi3 - ^ + 5^ v"'' 
X 7"'^^\ki,k2){k2 + 2{N-l)a,ki+2a\+h.c., 

ki ,k2 

which becomes the white noise in the thermodynamical 
limit because the off diagonal elements rapidly vanish 
with N. It can be checked that for larger separation one 
always gets the white noise for the same reason, i.e., the 
off-diagonal elements vanish with N. 

We finally calculate the negativity between nearest 
neighbour spins. Since their state is a mixture of a max- 
imally entangled state with identity, the overall eigenval- 
ues will be the same mixtures of the eigenvalues of the 
identity and the maximally entangled state. This is also 
true for the partially transposed state. It is therefore easy 
to see that the negative eigenvalues of the partially trans- 
posed state all have the absolute value of {S — l)/25^. 
The number of negative eigenvalues of the partially trans- 
posed total state is calculated to be {S^ — S)/2. The total 
negativity, which measures the amount of entanglement 
in the state [||, is hence S{S - 1)2/45^. We see that for 
large spin value, the total negativity grows linearly with 
the size of spin as claimed earlier. 

Proof of 2 and 3. The density matrix Po{Pe) describing 
the subset of all the odd (even) sites has the following 
form 

^-Nab 

Po = '2^ |fci,A:2,...,fcAr) 
fci ,k2,...,k]<j 

{kN + 2a, fci + 2a, . . . , kN^i + 2a\ + h.c. (7) 
Pe = 2-t--^ X! \ki+a,k2 + a,...,kN + a) 

ki ,k2- ..kN 

{k2,k3, . . . , kN, ki \ + h.c. 

Both of them become the white noise in the thermody- 
namic limit, which implies that there is no entanglement 



between any subset of odd (even) sites. However, the set 
of all the odd sites is maximally entangled with the set 
of all even sites. 

The entanglement between the subsets of odd and even 
sites can be seen already at the level of four sites. For 
instance, the state of the first four sites is given by 

P1234 = \{ab)i2, {ab) 34) {{ab) 12, (06)34! 

+ ^^\iab)23){{ab)23\®^ (8) 

+0{S^-'')\{ab)i3, (06)24) ((a6)i4, (06)23! + h.c. 

It is clear that the the sites 1 and 3 treated as one sub- 
system are entangled to the subsystem consisting of the 
sites 2 and 4 but there is no genuine multipartite entan- 
glement. The above formula can be easily generalized for 
an arbitrary subsets containing even and odd sites and it 
can be seen that as long as the size of the subsets is fixed 
the subset is not genuinely multipartite entangled in the 
thermodynamical limit. 

There are situations in which different dimer coverings 
have the same energy. This happens for example in the 
Majumdar-Ghosh Hamiltonian where we have near- 
est and next nearest neighbour interactions. Unless there 
is some broken symmetry mechanism each different cov- 
ering will contribute to the overall state with equal weight 
in the thermodynamic equilibrium. This means that is 
more appropriate to consider mixture of dimer coverings 
rather than their superpositions. 

We now consider mixtures a of this type in any spatial 
dimension 



where \ck) is fc-th dimer covering and M is the num- 
ber of all possible dimer coverings. For simplicity we 
present our result in two dimensions. The generalization 
to higher dimensions is straightforward. 

We consider an infinite two dimensional square lattice 
that is an union of two sub-lattices La and Lb. ^ site be- 
longing to the sub-lattices La have neighbours belonging 
to the sub-lattices Lb (coordination number is 4). The 
spin liquid state on the lattice Q is defined as the su- 
perposition of all possible dimer coverings between the 
sub-lattices La and Lb. The question arises what the 
ratios with which we should mix the coverings are. 

To answer this question let us imagine a situation when 
two arbitrary neighbours, say i and j , are in a maximally 
entangled state, i.e., they form a dimer. The rest of the 
sites can be covered by dimers and there are S of such 
coverings. For an infinite lattice S is, of course, infinite 
but we can assume that the lattice is very large in which 
case S is also large but finite. Suppose now that the 
same neighbours are not in a maximally entangled state, 
i.e., they form dimers with their other neighbours. As 
is illustrated in Fig. ([1]), it is easy to see that in each 



FIG. 1: This figure illustrates the mixture of four possible 
dimer coverings, {\ck),k = l,---,4}, for an infinite lattice 
site. The top box depicts one of the dimer coverings when 
the neighbors A and B are maximally entangled. The other 
three show some of the possible coverings when the state of A 
and B is maximally mixed (therefore disentangled.) We argue 
that the ratio of the number of the coverings when A and B 
is maximally entangled to the number of coverings when A 
and B is maximally mixed is 1 : 3 in the thermodynamic limit 
(just like in the figure). 

such case (there are three of them) the number of the re- 
maining dimer coverings S' is equal and approximately 
the same as in the previous case where the two neigh- 
bours i and j were in a maximally entangled state, i.e., 
T, K T,' . Therefore, for an infinite lattice the ratio of the 
dimer coverings including the sites i and j to the dimer 
coverings not including the sites i and j equals ^. The 
consequence of this is that after tracing out all the other 
sites the density matrix describing the neighbouring sites 
i and j {i and j can be, in fact, any neighbours) is given 
by a Werner state containing | of the white noise. Thus, 
for a system consisting of spin i there is no entanglement 
between the neighbours. The situation changes for larger 
spins S > 3 (recall that S — 2s + 1, where s is the spin) 
because of the higher amount of the white noise that 
can be admixed to a maximally entangled state without 
destroying entanglement Q . 

The similar reasoning can be applied to a "honey- 
comb" two dimensional lattice (coordination number 3). 
In this case each spin has only three neighbours which 
gives us the amount of the white noise in the Werner 
state between the nearest neighbours equal to |. For 
spin i we do not have entanglement (the state is on the 
verge of being entangled) whereas entanglement exists 
for larger spins, i.e, for S > 2. 

In higher than two dimensions we can apply exactly the 
same logic. All we need to do is to calculate the ratio of 
the number of coverings containing a dimer between two 
neighbouring sites and those that do not. This ratio is 
always a fraction , where i? is a finite number that is the 
function of the coordination number (for a simple regular 
lattice it equals to the twice of the spatial dimension). 
Therefore, in any dimension and any lattice structure 



FIG. 2: This figure depicts the arrangement of dimers in the 
BCS ground state in the momentum space. The pairs of 
points represent pairs of spin entangled electrons with mo- 
menta of the same magnitude but of opposite directions. As 
we see this arrangement is not isotrpoic, which is why our con- 
siderations no longer apply. In this state we have fermionic 
pairs condensation and yet each of the fermions is spin i. 
This is contrast to the isotrpoic case when the minimal spin 
required for entanglement in three dimensions is | [See text 
for detailed explanation]. 



entanglement will always exist for spins with magnitude 
higher than R; likewise it vanishes below this value. 

A possible test of our prediction that there is a crit- 
ical value of the spin below which entanglement does 
not exist could go as follows. Recently, condensates of 
fermionic atoms were observed in the laboratory [Sj]. In 
order to achieve condensation fermions have to form en- 
tangled pairs (so called Cooper pairs), which then behave 
like quasi-bosons. Entangled fermionic pairs are formed 
through scattering process the strength of which can be 
controlled experimentally by Feshbach resonance. Our 
calculation imply that only atoms with sufficiently high 
spins can form entangled pairs. In three dimensions we 
require atoms to have spin higher than |. The experi- 
ment by Regal et. at [8] uses potassium atoms, ^"X, 
whose spin is |. Were they to use atoms with spin of 
I or lower, we conjecture that no fermionic condensa- 
tion would be possible. Note that allowing long range 
dimers can only increase the spin value necessary for en- 
tanglement between any two sites on the lattice. This 
is because the ratio of the dimer coverings contributing 
to maximally entangled state between any two points on 
the lattice to the rest of the coverings decreases. 

It is important to emphasize that our assumptions do 
not need to hold in practice. For example, the BCS model 
of fermionic condensation [9| has a ground state where 
momenta of spins are coupled in opposite directions as 
is schematically shown in Fig. ([2]). Namely, spin dimers 
now exist only between specific points on the lattice. This 
affects the ratio of mixing and allows for entanglement 
with low spin in higher dimensions. This model is highly 
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non-isotropic in the sense of lattice points not having the 
neighborhood. 

Dimer coverings present very simple entangled states. 
Yet they are physically relevant as they frequently occur 
in nature. Their entanglement structure is simple. Here 
we show that overall entanglement is fully characterized 
only in terms of nearest neighbour entanglements. In 
spite of its simplicity, we encountered one surprising fact. 
Namely, in higher than two dimensions the existence of 
entanglement is dependent on the size of the spin in each 
lattice site. Many questions remain open. For example, 



what happens when a degree of anisotropy is introduced 
into dimer coverings? Also, when does long range entan- 
glement arise if we allow long range dimers? We hope 
that our work stimulates a number of other interesting 
directions of research. 
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